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1 Introduction 

This manuscript is the written version of the material I attempted to cover 
at "Critical Dialogues in Cosmology". My position is that although BBN 
is inevitable in an expanding universe filled with radiation and baryons, at 
present there is a conflict between the abundances of the light nuclides pre- 
dicted by standard, big bang nucleosynthesis (SBBN) and those inferred from 
current observational data. This challenge to SBBN provides opportunities 
for astronomy/astrophysics, cosmology and particle physics. Perhaps there are 
unaccounted-for errors in passing from the observational data to the derived 
abundances (astronomy/astrophysics). Perhaps there are unaccounted-for sys- 
tematic corrections in translating the derived abundances from "here and now" 
to "there and then" (astrophysics/cosmology). Or perhaps the tension between 
theory and observations is offering us a hint of new physics beyond the standard 
model (cosmology/particle physics). 

1.1 The Basics 

Which nuclei may be synthesized during BBN, and in what relative abundances, 
depends on the interplay among the nucleon density, the radiation temperature 
and the early expansion rate. For example, if the nucleon density were very 
small and/or the universe were expanding very rapidly when thermal energies 
were comparable to nuclear binding energies, few complex nuclei would have 
emerged from the big bang. In SBBN (isotropy, homogeneity, three flavors 
of light neutrinos, etc.) the yields depend on only one "free" parameter, 77, 
the ratio (at present) of nucleons to photons (rj — tib/Vi 7 ; r]w = 10 10 ?7). In 
Figure 1 are shown the predicted SBBN yields (Yp is the 4 He mass fraction and 
2/2) 2/3, 2/7 are the ratios by number to H of D, 3 He and 7 Li) as a function of 77. 
The present contribution of baryons to the total density, in units of the critical 
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density, is directly proportional to 77. For a present cosmic background radiation 
temperature of 7c br = 2.73K, 

n B = 0.0146^0/1^, (1) 

where the Hubble constant is Ho = 50/iso km/sec/Mpc. As may be seen in 
Figure 1, the relative abundances vary considerably with 77 so that SBBN is 
an overdetermined system in the sense that fixing one of the primordial abun- 
dances fixes i] and leads to predictions of the abundances of the three other 
nuclides. SBBN is an eminently testable theory! Consistency demands that the 
value of 7/ determined by, for example, the primordial abundance of D lead to 
predictions for the primordial abundances of 3 He, 4 He and 7 Li consistent with 
the observational data. 

1.2 Status Quo Ante 

Until recently, the program of confronting the predictions of SBBN with the 
observational data had led to increasing confidence in the consistency of the 
standard model (e.g., WSSOK For example, solar system and interstellar 
(ISM) observations of D and 3 He have been used in concert with Galactic evolu- 
tion models [|[ |[ f|, || (or, nearly model- independent "inventories" [0, @, ||) 
to bound the primordial D abundance from above, leading to a lower bound 
to 77 (7710 > 2.8) 0]. At the same time, observations of 4 He in low-metallicity, 
extragalactic HII regions ||, |l0[ [llj [l^] and of 7 Li in very metal-poor halo stars 
|H 0, |l| have led to upper bounds on their primordial abundances, provid- 
ing an upper bound to 77 (j]io < 4.0) [[I). This quantitative agreement between 
theory and observations for the abundances of the light nuclides, covering some 
nine orders of magnitude, provides strong support for the consistency of SBBN. 
Furthermore, it should be emphasized that the inferred value (range) of rj is 
in reasonable agreement with estimates of the baryon density derived from the 
dynamics of "luminous" (i.e., baryonic) matter (see, e.g., [T^ | and references 
therein) thus providing support for the extrapolation from the present universe 
to its earliest epochs. However, lest we fly too high on the sweet smell of suc- 
cess, it is sobering to recount the suggestions of problems with this quantitative 
comparison ]l7| ]. 

1.3 Hints Of A Crisis 

For SBBN the primordial abundances of D and 3 He decrease and that of 4 He 
increases with increasing nucleon-to-photon ratio (see Fig. 1), so that the lower 
bound to 77 inferred from observational upper bounds to D and 3 He leads to a 
predicted lower bound to the primordial 4 He mass fraction. Alternatively, the 
upper bound on Yp from the observational data leads to an upper bound to 77 and 
a lower bound to the primordial abundance of D (or D+ 3 He). Observations of D 
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and 3 He "prefer" a relatively large lower bound to T] while those of 4 He "favor" 
a relatively small upper bound to r\. This "tension" between D and 3 Hc on the 
one hand, and 4 He on the other, provides a hint of a problem (e.g., [[I], [| |l7|]). 
In Figure 2 this "crisis" is displayed. Superposed on the SBBN predictions for 
the abundances of D, 4 He and 7 Li are the 68% and 95% confidence level bands 
inferred from the observational data Jl7|] . Notice that the range of r\ delimited 
by D has (at 95%) no overlap with that inferred from 4 He. This is the "crisis" 
[|7j. Another way to visualize the crisis is provided by Figure 3 where the 68% 
and 95% CL ranges for r\ have been fixed by the inferred primordial abundances 
of D, 3 He and 7 Li, leading to predicted ranges for Yp which are to be compared 
to the upper bound derived from the HII region data fl2]| . Impressive as is the 
approximate consistency of SBBN, the standard model is seriously challenged. 

2 Three Possible Solutions 

The hint of a crisis outlined above bears some similarities to several crises which 
challenged the 19th century "standard model" (Newtonian mechanics/gravity) 
when solar system observations appeared to disagree with theoretical predic- 
tions. The "Uranus crisis" was real (accurate data) and was resolved in favor of 
the standard model leading to the discovery of something new: Neptune. The 
"Neptune crisis" was illusory (inaccurate data), the standard model remained 
consistent, and the discovery of Pluto was an "accident" . Most exciting, of 
course, was the resolution of the "Mercury crisis" . The crisis was real (the data 
accurate) and the standard model was superseded by General Relativity. Sim- 
ilar options appear available for the resolution of the current crisis. Perhaps 
the data used to derive the abundances of one or more of the light nuclides are 
tainted by unrecognized statistical or systematic errors (the Neptune crisis). Or 
perhaps the data are accurate but there are unforeseen astrophysical effects, 
unaccounted for in our extrapolation from the derived abundances to their in- 
ferred primordial values (the Uranus crisis). We would, however, be remiss to 
ignore the possibility (the Mercury crisis), that this conflict may be revealing 
evidence for new physics (cosmology and/or particle physics). Although the 
ultimate resolution of the current crisis may require some combination of the 
above suggestions, here I will consider them separately and identify possible 
examples of each. 

2.1 Inaccurate Data: 4 He ? 

Once the number of light neutrino flavors is fixed (N v = 3 for SBBN), the 
predicted 4 He mass fraction is a relatively insensitive function of r\ (see, e.g., 
Figs. 1-3). However, as 4 He is the second most abundant nuclide, Y may be 
determined very accurately in a variety of environments spread throughout the 
universe. One source of the current crisis is the relatively low upper bound 
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to r\ inferred from the primordial abundance (Yp — 0.232 ± 0.003 ± 0.005) de- 
rived from observations of recombination emission from hydrogen and helium 
in extragalactic HII regions |, Although there are large numbers of 
low metallicity HII regions which have been observed carefully, minimizing the 
extrapolation to the primordial abundance and leading to a very small statis- 
tical uncertainty in Yp (0.003), perhaps there are large systematic corrections 
along the path from the observed equivalent widths to the derived abundances 
which exceed the estimate of 0.005 For example, since neutral helium is 

unobservable, perhaps there is hidden neutral helium in regions where hydrogen 
is fully ionized. If unaccounted for, this would bias the results to abundances 
which are systematically low. The observers (clever as they are) have antic- 
ipated this possibility and have restricted attention to "high excitation" HII 
regions where observations of other ions and comparisons with models suggest 
this ionization correction is negligible. Indeed, for the low metallicity HII re- 
gions observed, the hard radiation spectrum from the very hot stars may even 
reverse this effect (neutral H where He is fully ionized) . Such an effect might 
even correlate with metallicity (downward correction to Y for low metallicity, 
upward correction for higher values) , leading to a reduction in the derived value 
of Yp . Collisional excitation (especially from the metastable level in Hel) could 
enhance the helium emission leading to an overestimate of Y. Most estimates 
suggest this effect is small, and the observers do, in one of several ways, try to 
correct for it. Perhaps they have overcorrected, leading to an abundance which 
is systematically too low. Again, since collisional excitation is temperature- 
dependent and metal-poor HII regions might be hotter, this systematic effect 
might correlate with metallicity. 

The bottom line here is that if Yp = 0.246 (rather than 0.232; see, e.g., 
[|l8|), the crisis would be resolved (in favor of a relatively higher value of 77 and 
relatively lower values of primordial D and 3 He and a higher value for primordial 
'Li). 



2.2 Uncertain Extrapolation: D ? 

Until recently, deuterium had only been observed "here and now" in the solar 



system |9[ and the (very local) ISM |2fj, |2lj . The older Copernicus UV data 
and the newer HST data lead to a reasonably accurate value for the abundance of 
ISM D (D/H = 1.6±0.2 x 10~ 5 ). Solar system data (meteoritic 3 He abundances, 
deuterated molecules in the giant planets, etc.) lead to a slightly larger (and 
slightly less accurate) abundance p9| . So far, so good. Although these data are 
local (in space and in time) , they are still of great value since the evolution of 
deuterium is simple. Other than the big bang, there are no astrophysical sites for 
the production of significant abundances of D [^2| , and D is destroyed whenever 
gas forms stars. Thus, yi (actually X2, the deuterium mass fraction) has only 
decreased since primordial nucleosynthesis ended (7/2P > 2/21SM ; actually, X2P > 
A^ism)- As a result, observations of deuterium anywhere, anytime provide a 
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lower bound to its primordial abundance and, therefore, an upper bound to rj. 

The bad news is that both the present ISM and the presolar nebula are 
evolved; they contain gas which has been processed through stars (where D is 
destroyed). To bound r\ from below requires that we bound j/2P (^2p) from 
above and this requires knowledge of Galactic chemical evolution. For a wide 
assortment of independent (but similar) chemical evolution models ||, ||, || 
p3| , designed to account for the observed age-metallicity relation, the observed 
gas/stars ratio, various abundance ratios (e.g., secondary to primary nuclei), 
etc., it is predicted that D is destroyed by a factor of 2 - 3 0. Given the strong 
dependence of j/2P on r\ (see, e.g., Figs. 1 & 2), such a modest destruction 
factor coupled with the accurate ISM abundance leads to a reasonably narrow 
range for r\ which, however, corresponds to a SBBN predicted abundance for 
4 He larger than that inferred from the HII region data (the "crisis"). Although 
reasonable, these upper bounds to primordial D are model-dependent. Could 
"designer" models for Galactic evolution be found which, while maintaining con- 
sistency with the wealth of observational data, nonetheless permit much larger 
D destruction? Without detailed models which have actually been confronted 
with the data, it is difficult to answer this question. Nonetheless, there are 
reasons to believe it will not be easy to find such models based on several nearly 
model-independent approaches which have been explored previously. 

2.2.1 The D+ 3 He Inventory 

One such approach to bounding D destruction has been to exploit the fact that 
when D is incorporated into stars and burned, it is first burned to 3 He p5| . 
The more resilient 3 He burns at a higher temperature than D, and so, while 
all the D is destroyed, some 3 He survives stellar processing. Indeed, in stars of 
all masses there are interior zones where hydrogen burning results in the pro- 
duction of new 3 He. Thus, in general, the more gas cycled through stars (and 
the more D destroyed), the more 3 He mig ht be expected §, §, [|| Unfor- 
tunately, the evolution of 3 He is very complex involving a competition between 
new production and the destruction and survival of the prestellar D+ 3 He. New 
production of 3 He is uncertain. If it is ignored current observations (ISM and/or 
solar system) of D and 3 He [|l9], |2(], ^l) may be used to provide upper bounds 
on primordial D and 3 He, nearly independent of the details of Galactic chemical 
evolution jl], ^, |t], ^ . All the unknown model-dependence of stellar and Galac- 
tic chemical evolution is contained in (53), the average 3 He "survival fraction". 
Such an approach was pioneered by Yang et al. Q and has been refined and up- 
dated by Steigman & Tosi and Hata et al. Q . For (53) > 1/4 and employing 
solar system D and 3 He and ISM D data, Hata et al. |§] derive for SBBN the 
best fit values (95% CL): (D/H) P = 3.5±?;| x 1CT 5 and r] 10 = 5.0±H. This 95% 
CL range for 77, derived from D and 3 He alone, is considerably higher than the 
previously preferred range (WSSOK) , exacerbating the tension between D and 
4 He since it predicts (for SBBN) Yp — 0.247 ±0.004, far in excess of the value in- 
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ferred from the HII region observations [jT2[ . This strong disagreement could be 
ameliorated if (53) is smaller than the value (1/4) usually adopted [jl], ^, |?|, |L p?J . 
For (53} < 0.1, considerably less 3 He survives, permitting more destruction of 
D, consistent with a larger primordial abundance (corresponding to lower values 
of 77 and Yp). Production of new 3 He, if any, would have an effect similar to 
that of increasing (33), restricting primordial D and 3 He to lower abundances 
leading to higher values of 77. There are reasons to suspect that current stellar 
models may overestimate the production of 3 He || |8[ |9|, [k| [H). But, even if 
total destruction in stars less massive than 2.5M Q is assumed, (173} > 0.3 for gas 
which has been cycled through one (and only one) generation of stars [|| . Thus 
to reduce (33) would seem to require a model where gas was efficiently cycled 
through several generations of stars. As outlined in the next section, there may 
be problems with such models. 

2.2.2 Constraints On D Depletion From Metallicity And Gas/Stars 

As mentioned above, to relax the upper bound on the primordial abundance of D 
(and, correspondingly, the lower bounds on 77 and Yp) seems to require chemical 
evolution models which are efficient in cycling gas through stars. 3 Hc provides 
one possible constraint on such models which, however, may be plagued by un- 
certainties in stellar modelling. The observed metallicity provides another such 
constraint ||^, p^| . The more gas cycled through stars, the higher the metallic- 
ity in young stars and newly returned gas. It will be a challenge for "designer" 
models of Galactic evolution to accomodate large D destruction while avoiding 
overproduction of the heavy elements. Perhaps this might be accomplished by 
expelling (via winds or superbubbles?) the metal-enriched gas while retaining 
the D-depleted gas; it remains to be seen whether realistic models of this type 
can be found. 

Another challenge to such models is the observed ratio of mass in gas to that 
in stars (e.g., |3J]). Is the relatively high observed gas/stars ratio consistent with 
models which efficiently cycle the gas through stars (e.g., see |||])? Such models 
might require infall to replenish the gas supply, but if the infalling gas consists 
largely of unprocessed material, the ISM D abundance is driven back towards 
its primordial value. 

Nonetheless, given our ignorance of Galactic chemical evolution, we cannot 
dismiss the possibility that D has been destroyed by a large factor (~ 5 - 10) 
between the big bang and the present epoch. If so, lower values of 77 might be 
compatible with ISM and solar system observations of D (and 3 He), leading to 
lower predicted values for Yp , consistent with the HII region data. Perhaps the 
astrophysics is at fault. Recent data from high rcdshift (z), low metallicity (Z) 
QSO absorbing clouds present ambiguous support for this possibility (see Sec. 
3). 
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2.3 New Physics: A Massive Tau Neutrino ? 



Due to the gap at mass-5 and the strong binding of 4 He, once BBN begins most 
available neutrons are burned to helium-4. As a result, although Yp is relatively 
insensitive to 77, it is closely tied to the neutron-to-proton ratio at BBN. In 
turn, n n /n p is regulated by the competition between the weak interaction rates 
(n <-> p) and the universal expansion rate (H). H is fixed by the total energy 
density at BBN and, for SBBN, is dominated by the contribution from the 
extremely relativistic particles present (7, e , v e , v^, v T ). In units of the photon 
density, for SBBN, PtcttVay = 43/8. If, due to "new physics", H is changed 
from its SBBN value, n n /n p at BBN will be modified leading to a different 
Yp versus 77 relation. One such possibility is that the standard cosmology may 
be modified due to a variation in G, the Newtonian gravitational "constant" . 
Another is that the particle physics content of the early universe could change 
due to the presence of additional light neutrinos or light scalars, or if the tau 
neutrino were very massive (> ~ 10 MeV) and unstable p5| . It is convenient to 
parameterize such changes by N„ p6| , the "effective number of equivalent light 
neutrinos" (complementary to the number of standard model neutrinos probed 
by collider experiments) where N„ is defined by, 

Ptot/pIot N = 1 + 7(JV„ - 3)/43 (2) 

For N„ > 3, the early universe expands more quickly, leaving behind more 
neutrons to be incorporated into 4 He during BBN, and vice-versa. To a good 
first approximation, Alp as 0.0l(N v — 3). Since one aspect of the crisis is that 
SBBN consistent with the inferred primordial abundances of D, 3 He and 7 Li 
predicts a value for Yp which is larger than that derived from the extragalactic 
HII regions by AY P w 0.01, N v w 2 could reestablish consistency @ §]. Of 
the myriad possibilities, a massive (unstable) tau neutrino provides one option 
for reducing N v from its SBBN value of 3 to one closer to 2 [§5| . Current collider 
experiments are capable of exploring the interesting mass range (> ~ 10 MeV) 
and either ruling out this option or confirming it. 



3 High-z, Low-Z Deuterium 

As noted earlier, the relatively strong dependence of the SBBN predicted abun- 
dance of D on 77, coupled with the simple evolution of deuterium, identifies D/H 
as the ideal baryometer. It has been anticipated that observations of deuterium 
in high redshift (nearly primordial), low metallicity (nearly unevolved) QSO 
absorption systems would resolve the current confusion, relieving the tension 
between theory and observations. Unfortunately, the few such cases identified 
to date have added to the confusion, and rather than relieving the tension, have 
led to a rupture. As described eloquently by Hogan |37j], he favors the "high-D" 
results (y 2 p = 19 ± 4 x 10~ 5 ) §|, |§ while I'm sure Tytler fh], |IJ would argue 
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for the "low-D" values (y 2 p = 2.4 ± 0.3 ± 0.3 x 10~ 5 ). As may be seen from 
Figure 4 p2| , if either of these values is correct (rather than something between 
the two) then the "true" value of 77 is actually outside the old concordance range 
(2.8 < r/io < 4.0). For example, for "high-D", 1.3 < rj w < 2.7 (95%CL), while 
for "low-D", 5.1 < 7710 < 8.2 (95%CL) @. For "high-D" the SBBN crisis 
dissipates, since the predicted abundances of 4 He and 7 Li (Yp = 0.234 ± 0.002, 
2/7P = 1.5 ±0.6 x 10 -10 ) are in excellent agreement with their inferred primordial 
values |p3[ . In this case the challenge is in the chemical evolution court since 
the comparison of the high primordial D abundance with the low ISM value 
requires that deuterium should have been destroyed by a factor of ~ 13. In 
contrast, for "low-D", the SBBN crisis is exacerbated @ (Yp = 0.249 ±0.001, 
J/7P = 4.7 ± 0.7 x 10~ 10 ), while the D evolution (destruction by a factor of 
~ 1.6) is consistent with "normal" Galactic chemical evolution models (e.g., 
[|| ||, gj. As always (and as it should be) we await more data to resolve the 
current conundrum. 



4 A Different Path To The Baryon Density 

Given the current ambiguous state of affairs it may be of some value to explore 
non-BBN pathways to the baryon density. Here x-ray clusters may play a valu- 
able role. To the extent that rich clusters of galaxies provide a "fair sample" of 
the universal baryon fraction, the cluster baryon fraction may be used in con- 
cert with dynamical determinations of the total density (in clustered matter) to 
infer the baryon density. 

n B = f B n ; 77 10 = 273/sfi/i 2 (3) 

In equation (3), H = 100ft, = 70 ± 10 |Q and £1 is the density of clustered 
matter (= fie + Hcdm for "open" and "lambda" CDM models) in units of 
the critical density. Within the context of "open" and "lambda" CDM models 
observations of large scale clustering constrain the combination £lh « T « 0.25± 
0.05 [jt5| . From x-ray, optical and "mini-lensing" studies of rich clusters, /g 1 « 
/hg + /gal w (0.07 ± 0.01)ft- 3 / 2 + (0.02 ± 0.01) so that [§§ 7710 « 6.7 ± 1.6. 
Actually, further bounds on f2 and Hq from large scale velocity flows |^7j and 
the age of the Universe tend to increase f2 and reduce Hq from the rough 
estimate presented here, leading to a somewhat higher range for r\ |^6j . This is 
true also for "mixed" hot plus cold dark matter models where Q > 1 — fin dm > 
0.7; even for Hq = 55 ± 10 [Q, 7710 > 11 ± 2. Thus, in the context of a variety 
of CDM models, the preliminary results of this "dynamical" approach to the 
baryon density suggest a high value for rj, consistent with that inferred from 
solar system and ISM D and 3 He ||], favoring the low-D/high- 4 He option. 
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Figure Captions 

Figure 1. SBBN predicted primordial abundances of 4 He (mass fraction 
F P ), D (y 2P = (D/H) P ), 3 He (y 3P = ( 3 He/H) P ) and 7 Li (y 7P = ( 7 Li/H) P ) as a 
function of 77, the nucleon-to-photon ratio. This graph has been provided by D. 
Thomas. 

Figure 2. SBBN predictions (solid lines) for Yp, y 2 p and y?p along with 
their theoretical uncertainties (la) estimated via Monte Carlos (dashed lines) 
| ji"7[ . Also shown are the regions constrained by the observations at 68% and 
95% (shaded regions and dotted lines respectively). 

Figure 3. Theoretical (SBBN) predictions for the primordial hclium-4 mass 
fraction (1th) at 68% (shaded region) and 95% (dotted curve) for i] constrained 
by the inferred primordial abundances of D, 3 He and 7 Li. Also shown are 
the 68% and 95% bounds to Yp inferred from observations of low-metallicity, 
extragalactic HII regions (Yobs). The absence of overlap between Y ODS and Y t h 
is the "crisis" . 

Figure 4. As for Figure 2 but, with the deuterium abundance inferred from 
the two sets of QSO absorption data Jf2"| . 
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